Droplets in the two-dimensional ±J spin glass: evidence for (non-) universality 
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Using mappings to computer-science problems and by applying sophisticated algorithms, one 
can study numerically many problems much better compared to applying standard approaches like 
Monte Carlo simulations. Here, using calculations of ground states of suitable perturbed systems, 
droplets are obtained in two-dimensional ± J spin glasses, which are in the focus of a currently very 
lifely debate. Since a sophisticated matching algorithm is applied here, exact ground states of large 
systems up to L 2 — 256 2 spins can be generated. Furthermore, no equilibration or extrapolation to 
T = is necessary. Three different ±J models are studied here: a) with open boundary conditions, 
b) with fixed boundary conditions and c) a diluted system where a fraction p = 0.125 of all bonds 
is zero. For large systems, the droplet energy shows for all three models a power-law behavior 
En ~ L B ° with #fj < 0. This is different from previous studies of domain walls, where a convergence 
to a constant non-zero value (#dw = 0) has been found for such models. After correcting for the 
non-compactness of the droplets, the results are likely to be compatible with fo ~ —0.29 for all 
three models. This is in accordance with the Gaussian system where #d = —0.287(4) (i> ~ 3.5 
via v = — 1/#d). Nevertheless, the disorder- averaged spin-spin correlation exponent r\ is determined 
here via the probability to have a non-zero-energy droplet, and rj « 0.22 is found for all three models, 
this being in contrast to the behavior of the model with Gaussian interactions, where exactly r\ = 0. 

PACS numbers: 75.50.Lk, 02.60.Pn, 75.40.Mg, 75.10.Nr 



Monte Carlo simulations [j], Q and related approaches 
are common ways to study physical systems numerically, 
in particular for the usual case where no exact analytic 
solution can be provided. So far, it has only occasionally 
been recognized that by using mappings to computer- 
science problems and by applying sophisticated algo- 
rithms one can obtain in many cases results which are 
superior in comparison to using standard algorithms. 

Here, Ising spin glasses are considered, which are the 
most-frequently studied systems in statistical physics 
@) 0, H, H • However, despite more than two decades 
of intensive research, many properties of spin glasses, 
especially in finite dimensions, are still not well un- 
derstood. For two-dimensional spin glasses it is now 
widely accepted that no ordered phase for finite tem- 
peratures exists @, H, H, E3] m this case. Nevertheless, 
the ±J model with a bimodal distribution of the inter- 
actions is in the center of a currently very lively debate 
H, EH, El EI EI EE El, El El, EH m EH , in particular 
whether the behavior is equivalent to that of the model 
with a Gaussian distribution of the interactions. One 
central question is whether the correlation length, when 
approaching T c = 0, diverges algebraically [22J ~ T~ v ', 
as for the Gaussian model, or like an exponential [23} 
r~> T~ 2 e~ c / T , formally equivalent to v = 00. In partic- 
ular, Ref. El claims evidence for a power-law divergence 
with the same critical exponent v w 3.5 as the Gaus- 
sian system, and that the spin-spin correlation exponent 
i] = 0. This exponent is defined via [(SiSi + i) 2 ]j ~ l~ v , 
[...}, j and (...) being the averages over the quenched dis- 
order and the thermal average, respectively. 

Unfortunately, most of the above cited work is based 
on finite-temperature calculations, in particular Monte 
Carlo simulations, hence an extrapolation T — ► T c = is 



necessary. Also the systems are restricted in most cases 
to rather small sizes L < 64. Only in the case, where par- 
allel tempering Monte Carlo simulations 16j, [21 1 or the 
worm algorithm [IE] have been used, sizes L < 128 could 



be considered. Consequently, in Ref. [21| it has been shown 
that the presently available finite-temperature data does 
not allow to draw final conclusions. In Ref. LLJ also ex- 
act T — properties are calculated (up to L — 128 for 
Lx L systems and up to L = 64 for 9L x L systems), but 
only the spin-spin correlation could be obtained by this 
approach. 

Here, we go much beyond t he p revious work. Exact 



ground-state (GS) calculations 24 1 are applied, which al 



27j for large systems like 



lows to obtain GSs [H, [H, [H, 
L = 480. The method relies on mapping the GS calcu- 
lation to a graph-theoretical problem and using sophis- 
ticated algorithms developed in computer science. By 
using suitable perturbations of the original systems, one 
can go beyond pure GS calculations and study excitations 
like domain walls (DWs) [H, HU and droplets [H [27[. 
Hence, large systems in exact equilibrium can be investi- 
gated and no extrapolation to T = is necessary. 

In Refs. [lfj, l25j, l26j, |27| this approach has been used to 
show that for the two-dimensional spin glass with Gaus- 
sian disorder all assumptions made by the droplet theory 
[28L I29I l3(i 31 j are fulfilled. In particular the energy 
scaling of the basic excitations DWs and droplets fol- 
lows power laws E dw ~ L 6d ™ and E D ~ L e °, respec- 
tively, with the same universal value 9 = 6^ = 9d w , 
which is related [32| to the correlation-length exponent 
via 6 = -1/v = -0.287(4). On the other hand, for the 
± J model, the average energyof DWs approaches a con- 
stant for large system sizes [10( L > 128, i.e. #dw = 0. 
This appears to be, via 9 — compatible with an 
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exponential divergence of the correlation length, if either 
the low-temperature behavior is dominated by DWs, or if 
the scaling of the energy of droplets, which are expected 
to be the dominant excitations, is the same as for DWs. 

In this work, droplet excitations for three different 
models with discrete distributions of the interactions are 
calculated using an approach which is based on exact GS 
calculations. This approach allows to consider large sys- 
tems up to L — 256, no extrapolation of the temperature 
and no equilibration are necessary. The main result is 
that the scaling behavior is different from the scaling of 
DWs, i.e. 9o < 0. This is compatible with a power-law 
divergence of the correlation length, as recently claimed 
[ill. After correcting for the non-compactness of the 
droplets [33| , the results seem to be in all three cases the 
same as for the Gaussian model, hence Gaussian and dis- 
crete models appear to be in the same universality class 
regarding the low-temperature behavior. Furthermore, 
in this work the scaling of the spin-spin correlation func- 
tion is studied, resulting in r\ w 0.22 for all three models. 

The Hamiltonian which is studied here is the usual 
Ising spin glass model: 



n 



(1) 



where the spins Si = ±1 lie on the sites of a square lattice 
with N = L 2 sites, the bonds Jij couple nearest-neighbor 
sites on the lattice. Jy = ± J with equal probability for 
the quenched realizations of the disorder. Here systems 
with either open boundary conditions (be) in all direc- 
tions, for the initial GS calculation, or fixed be always are 
studied. For the latter case, also diluted samples are con- 
sidered, where each bond is set to zero with probability 
p = 0.125; 

Here, droplets are considered, as introduced by 
Kawashima [33j , which are lowest-energy excitations with 
respect to the GS. They consist of a connected cluster of 
spins, which include a certain pre-selected spin, here a 
center spin of the system. The spins at the boundary are 
fixed to their GS orientations. The energy scaling of these 
droplets is expected to follow a L e ° scaling, where (33| 
9'y) = 9^>Dy jD is related to the volume fractal dimension 
Dy of the droplets, the dimension D = 2 of the system 
and the usual droplet exponent 9r>- The main approach 
used here is based on mapping the GS calcu- 

lation to the minimum-weight perfect matching problem 
and using sophisticated matching algorithms from graph 
theory. For details, please see the pedagogical descrip- 



tion in Ref. |24|. The droplet calculation of each disorder 
realization is based on a sequence of 2L suitable modifi- 
cations of the disorder, each time followed by a GS cal- 



culation. The details of the algorithm [37| are described 
extensively in Refs. IS S3]- 

Minimum-energy droplets have been obtained using 
the aforementioned approach for system sizes L = 6 to 
L = 256 (L = 160 for the fixed be). All results are 



averages over many disorder realizations, the number of 
realizations is between 20000 for small sizes and 5000 for 
the largest sizes. 
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FIG. 1: Average droplet energy Eo as a function of system 
size in a double logarithmic plot. The solid line sows a fit to 
a power law with correction to scaling, while the dashed lines 
show the results of fits to a simple power laws. 

In Fig. Q] the average droplet energy Ed is shown as 
a function of system size L. For large system sizes, a 
power-law behavior is visible. This corresponds 32j to 
a power-law divergence of the correlation length as ob- 
served recently by Jorg et al [la ], as in contrast to the 
behavior of DWs for the ± J model and other models ex- 
hibiting a quantized energy spectrum. A fit to the func- 
tion AL 9 ' D for L > 32 yields 6£, = -0.244(6) for open be 
with a good quality of the fit [H,[3§| Q = 0.32. Similarly, 
9' D = -0.250(6) (Q = 0.82) are obtained for fixed be and 
9' D = -0.295(6) (Q = 0.64) for the diluted systems. 

To obtain the droplet exponent 9, the geometric prop- 
erties of the droplet volumes are considered next. In 
Fig. [2] the volume V of the droplets is shown as a func- 
tion of the system size L for the three different mod- 
els. Note that due to the degeneracy of the models, 
droplets with many different values for the volume are 
possible. Unfortunately, unless through complete enu- 
meration, no algorithm to sample GSs or droplets with 
the same weight/probability is known. Here, the degen- 
eracy is broken, by selecting the droplets with the small- 
est surface [4ol | , hence the behavior of the volume is not 
controllable and can give only a rough idea of the true be- 
havior. When fitting a power law ~ L Dv to V(L) of the 
open be model, one obtains Dy = 1.81(1), Dy — 1.94 for 
fixed be, while the behavior of the diluted system is com- 
patible with Dy = 2. Via considering 6>d = 9' D D/Dy, it 
appears likely that 6>d ~ —0.29 universally for all three 
models. Since the low-temperature behavior is domi- 
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FIG. 2: Volume V of the droplets as a function of the system 
size, respectively. Error bars are smaller than symbol sizes. 
The lines show the results of power-law fits of ~ L Dv to V(L) 
with Dv = 1.81, Dy = 1-94 and D v = 2 for open be, fixed 
be, and the diluted system, respectively. The data for the 
diluted system has been shifted up by a factor of 2, for better 
visibility. 



nated by droplets, and not by DWs, this would mean 
the ± J models not only show a power-law behavior for 
the divergence of the correlation length like the Gaussian 
system, but exhibit the same value v = — 1/#d w 3.5 for 
the correlation-length exponent. In any case, even if v 
could still be slight different for Gaussian and ± J Ising 
spin glasses, the results show unambiguously that the 
scaling behavior of droplets and DWs is different in the 
±J case, the main reason 12J being the discreteness of 
the spectrum of excitations. This i s oppo sed to the as- 



j oppo sed t 

sumptions of the droplet theory |28j, 2iJ 53, 3l| , were the 
behavior of all types of excitations is governed by one 
single exponent 9 — 6*d = Odw 

Finally, to determine 77, here the relation [29[ 
[(SiSi + i) 2 ]j = p(L) is used, where p(L) is the probability 
to have a droplet with non-zero energy. In Fig. [3]p(X) 
is shown as a function of system size. According droplet 
theory [29} one expects p(L) ~ L~ r> . When fitting for 
system sizes L > 32, we obtain in all three cases 77 close 
to = 0.22 with good qualities of the fit, showing the uni- 
versality of this result. This makes it likely that also the 
behavior of the average droplet energy is universal, hence 
the differences observed in Fig. Q] are indeed through dif- 
ferent fractal dimensions Dy- Note that if one related 
77 to the probability that a non-zero-energy cross-system 
DW exists, one would assume 77 = since the DW en- 
ergy settles at a constant value for large systems ficj . 
i.e. #dw = 0. Nevertheless, two spins are uncorrelated if 
there exists any zero-energy excitation separating them, 



FIG. 3: Probability p(L) that a droplet has non-zero energy 
as a function of system size in a double logarithmic plot, for 
open be, fixed be and the diluted system, respectively. The 
lines show simple power laws ~ L~ v , with 77 = 0.22 in all 
three cases. 



hence one has to consider droplets as well, as performed 
here. 



The result is compatible with r\ = 0.21 obtained [41 1 
from T = transfer-matrix calculations of the correla- 
tion function for L < 12. On the other hand, the claim 
7] = of Ref. [HI, which is based on an extrapolation 
T — > for small systems L < 64, x2 is clearly ruled out. 
Even more, fits of the actual data obtained in Ref. [HI 
yielded, depending on the fits, also values 77 > 0. In Ref. 
16l a value 77 = 0.138 was found via Monte-Carlo sim- 
ulations at finite but low temperature for system sizes 
L < 128, but it was mentioned by the authors that their 
results are compatible with a "large range of 77 values" . 
Interestingly, at higher temperatures, a higher effective 
exponent rj e ff w 0.2 was observed by the authors. A value 
7/ = 0.14(1) has been obtained in Ref. [lj, but that result 
is obtained from studying only four different system sizes, 
and the result depends on the assumption that the corre- 
lation length at T — diverges as L 3 ^ 2 . Note that for the 
Gaussian model at exactly T = 0, due to the uniqueness 
of the GS, we have 77 = 0, hence regarding this quantity, 
the two classes of models look non-universal. 

One can understand by a simple scaling argument, 
why DWs and droplets can behave differently in the ± J 
model. We denote by qi the probability that a zero- 
energy DW exists in a system of size I. Since this 
probability approaches [10( a finite value q for about 
L > 100, we can assume qi ~ q for simplicity. We look at 
the system at different scales I — l®, Iq, Iq, . . . , Z§ = L 
(Iq > 1 arbitrary, k — lnL/ln/o), where we can as- 
sume the different DWs are independent. A non-zero 
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droplet exists only if on all scales no zero-energy closed 

DW exists, i.e. p(L) = \ \, r /; w q k = L~ v with 

— r\ = lng/ln^o < 0. Hence, the probability for non-zero 
energy droplets decreases with a power-law (hence the 
mean droplet energy), while the probability for non-zero 
energy DWs (and the mean DW energy) saturates for 
L — > co. 

To summarize, droplet excitations for three different 
variants of the two-dimensional ±J Ising spin glasses 
were studied. Here an advanced methodology from graph 
theory is used, based on mapping the GS calculation to 
the minimum- weight perfect matching problem, using so- 
phisticated matching algorithms from computer science 
and studying sequences of suitable modified realizations 
of the disorder. This allows to treat at T = 0, without 
need for an extrapolation of the temperature, large sys- 
tems up to L = 256 exactly. The average droplet energy 
shows a clear power-law behavior with exponent 8'^ < 0. 
It appears likely that when taking the non-compactness 
of the droplets into account, the same droplet exponent 
D k —0.29 emerges for all three models studied here, 
and hence the same value v = — 1/#d ~ 3.5 as for the 
model with Gaussian disorder is obtained. On the other 
hand, the value r) rs 0.22 for the exponent describing 
the decay of the spin-spin correlations is obtained for all 
three models, but this is clearly different from the Gaus- 
sian model, where 77 = 0. Remarkably, #d 7^ $dw = can 
be explained by a simple scaling argument. 
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